in 
o 



- 1—^ 

X 



UTAP-426 
astro-ph/0210502 

Supernova relic neutrinos and observational implications for 

neutrino oscillation 

Shin'ichiro Ando 1 '^ and Katsuhiko Sato 1 ' 2 

1 Department of Physics, School of Science, 
the University of Tokyo, 7-3-1 Hongo, 

Bunkyo-ku, Tokyo 113-0033, Japan 
Research Center for the Early Universe, 
School of Science, the University of Tokyo, 
7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 
(Dated: February 2, 2008) 



m 

O 2 

o ■ 

8 



CN ■ Abstract 

We investigate the flux of supernova relic neutrinos (SRN) for several neutrino oscillation models 



with parameters inferred from recent experimental results. In the calculation, we adopt the realistic 



' time- dependent supernova density profile, which is very different from the static progenitor profile 

O : 

owing to shock propagation. The Earth matter effect is also included appropriately using realistic 
density profile of the Earth. As a result, these two effects are found to induce the flux difference 

o 

by a few % in the detection energy range {E v > 19.3 MeV). We also set 90% C.L. upper limit on 
SRN flux for various oscillation models using the recently released observational result by Super- 



Kamiokande (SK). In the near future, further reduced upper limit is actually expected, when the 



current SK data are reanalyzed using some technique to reduce background events against SRN 
signals. It is expected that the reduced upper limit is sufficient to provide useful implications for 
neutrino oscillation. 
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I. INTRODUCTION 



A core-collapse supernova explosion produces a number of neutrinos and 99% of the 
gravitational energy is transformed to neutrinos. It is generally believed that the core- 
collapse supernova explosions have traced the star formation history in the universe and 
have emitted a great number of neutrinos, which should make a diffuse background. This 
supernova relic neutrino (SRN) background is one of the targets of the currently working 
large neutrino detectors such as Super-Kamiokande (SK). Comparing the predicted SRN flux 
with the observations provides potentially valuable information on the nature of neutrinos 
as well as the star formation history in the universe. This SRN background has been 
theoretically discussed in a number of previous papers 00 

On the other hand, there are observational constraints on the SRN flux. Most recently, 
SK collaboration set an upper bound of 1.2 v e cm -2 s _1 for the SRN flux in the energy region 
E v > 19.3 MeV [3j. (This limit can constrain only the flux of u e , since that flavor is the most 
easily detected by SK.) It is two orders of magnitude lower than the previous one obtained 
by Kamiokande II |4J , and is the same order as some typical theoretical predictions . For 
example, Ando et al. [3] (hereafter AST) predicted that the total SRN flux integrated over 
entire energy was 11 cm -2 s -1 , while the corresponding SK limit calculated with the AST 
spectral shape is 31 cm -2 s -1 . Since the theoretical calculations contain many ambiguities 
such as the supernova rate in the universe and neutrino spectrum from each supernova, 
this severe observational SRN limit can provide a number of valuable information on the 

n 

various fields of astrophysics and cosmology (e.g. Ref. |5|). Further, in the near future, it is 
expected that the upper limit will be much lower (about factor 3) when the current SK data 
of 1,496 days are reanalyzed using some technique to reduce background against detection. 
In that case, the SRN signal might be detected and the very severe constraint not only on 
the star formation history but also on the nature of neutrinos might be obtained. 

Thus, it is obviously important and very urgent to give more precise prediction for the 
SRN flux and event rate. For that reason, in this paper, we investigate the SRN flux 
using the most realistic models to date. In particular, we include the new features which 
have not been considered in all of the past studies Q, Q], illustrating them below. First, 
when we calculate the neutrino conversion probability in supernova, we adopt the realistic 
time- dependent density and Y e profiles, which are calculated by the Lawrence Livermore 
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group During the neutrino burst (~ 10 sec), the shock wave propagating the supernova 
matter changes density profile dramatically, and it is expected to affect the adiabaticity 
of resonance points [7]. Second, we consider the Earth matter effect. (For that effect on 

and references therein.) 



the future Galactic supernova neutrino burst, see Refs. 
Since SRN come from all the directions, half of them pass through the Earth matter, and 
it is also expected to change the SRN spectrum. Finally, neutrino oscillation with inverted 
mass hierarchy (777,3 <^ m i) is also investigated. In that case, the resonance also occurs in 
anti-neutrino sector, and it is expected that the SRN spectrum would be quite different from 
that in the case of normal mass hierarchy. Further, we repeat the discussions given in the 
recent SK paper 3] and obtain the 90% C.L. upper limit for the SRN flux estimated using 
various oscillation models. 

This paper is organized as follows. In Section m we give the models of original neutrino 
spectrum and supernova density profile. We also illustrate the neutrino oscillation models 
investigated in this paper, and qualitative behavior of flavor conversions for these models. 
In Section HTH we show the calculated SRN fluxes, and in Section HV| we compare the result 
with AST calculation. We also compare them to the observational upper limit from SK and 
discuss the validity of various oscillation models. Model dependence of our calculations are 
also discussed in the same section. 

II. ADOPTED MODELS AND NEUTRINO CONVERSIONS 
A. Supernova and Earth models 

As for the neutrino spectrum at production, we use that calculated by the Lawrence 



Livermore group 



llj . The mean energies are different between flavors, such as 

(E Pe ) ~ 16 MeV, (E„„) ~ 22 MeV. (1) 



This hierarchy of mean energies is explained as follows. Since 7> x 's interact with matter only 
through the neutral-current reactions in supernova, they are weakly coupled with matter 
compared to u^s. Thus the neutrino sphere of z/^'s is deeper in the core than that of z/ e 's, 
which leads to higher temperatures for u^s. Therefore, flavor conversions are expected to 
enhance the mean energy of z/ e 's. However, recent studies of neutrino flux and spectra 
formation in a supernova core (e.g. Ref. [12]) have shown that average v x energy exceeds 



the average v e energy by only a small amount, 10% being a typical number. We discuss this 
new aspect in Section HV CI 

We calculate neutrino conversion probabilities in supernova using the realistic time- 



dependent density and Y e profiles, w. 



rich are also calculated by the Lawrence Livermore 



group 



Among the past studies , AST first estimated the effects of neutrino oscilla- 
tion quantitativel y ho wever, they calculated conversionjprobability using static pre-collapse 
progenitor model |l3| . As the recent study indicates pj, the density profile changes dras- 
tically during neutrino burst (~ 10 sec) owing to shock propagation in supernova matter. 
When the shock propagates through the regions where matter-enhanced neutrino flavor 
conversion occurs (so called resonance point), it is expected to affect the adiabaticity of the 
resonance. (For the resonance and its adiabaticity, see the next subsection for details.) We 
show the pY e profile of supernova at 0.5, 2, 5, 10, and 15 seconds after bounce in the upper 
panel of Fig. As is seen clearly, the resonance takes place more than once at > 2 sec, and 
this is not expected when we use the static progenitor model. (As is discussed in the next 
subsection, the resonance for the anti-neutrino sector occurs only in the case of inverted 
mass hierarchy, or INV model in our notation.) 

Further, we include the Earth matter effect using realistic density profi 
This effect has not been considered in the past studies including AST 



Q. 



e of the Earth 

on the other 



hand, it has been discussed in the case of the future Galactic supernova neutrino burst by 
various authors (see Refs. 8j, l! llOl and references therein) . Although the effect is expected 

nnn 

to be very small for anti- neutrinos |8J, 121, LLUJ , we include it because the half of the SRN pass 
through the Earth matter. 



B. Neutrino conversions 



atest results which are updated recently by 



' r, l llfij - VVo adopt only large mixing angle 



As for oscillation parameters, we adopt the 
the SK and SNO solar neutrino observations 
(LMA) solution, whose parameters are 

Am 2 12 = 5.0 x 10~ 5 eV 2 , tan 2 12 = 0.42. (2) 

On the other hand, from the atmospheric neutrino experiments |l7], we use the values 

|Am 2 3 | = 2.8 x 10" 3 eV 2 , sin 2 26 23 = 1.0, (3) 



in our calculations. Whether Am 2 3 = mj — m 2 is positive (normal mass hierarchy) or 
negative (inverted mass hierarchy) is not known well, although the future Galactic neutrino 
burst is expected to provide useful information on the mass hierarchy (j| 0, [3| and actually 
several implications for it have been obtained with the neutrino burst from SN 1987A [19]. 
For # 13 , which is not also sufficiently constrained, we adopt two large and small values, 

sin 2 20i3 = 0.04, 1.0 x 10~ 6 , (4) 



both of which satisfy the upper bound from reactor experiment [2fJ. (Investigating neutrino 
oscillation using intermediate values for i3 is beyond the scope of this study, however, it is 
given in Ref. [6[ for the future Galactic supernova neutrino burst.) From these discussions, 
we adopt four parameter sets, named as NOR-S, NOR-L, INV-S, and INV-L, where NOR 
and INV represent the normal and inverted mass hierarchy respectively. The suffixes -L and 
-S attached to NOR and INV stand for large and small #13, respectively. In addition, we 
also investigate in case of no oscillation, for comparison. 

Now, we qualitatively discuss the behavior of anti-neutrino conversions during propaga- 
tion in supernova matter 1 (see e.g., Ref. for details). 

We consider, first, the case of normal mass hierarchy. The state of v e produced at 
deep in the core is coincide with mass eigenstate i>\, owing to large matter potential. This 
state can propagate to the supernova surface without being disturbed by the level crossing 
between different mass eigenstates (it is said that there are no resonance points). Thus, v e 
at production becomes V\ at the surface of the supernova and the observed v e flux is given 
by 



Fp e — \U e i\ F Pl + \U e2 \ Fp 2 + \U e3 \ Fp 3 

= \U el \ 2 Fl + (l-\U el \ 2 )Fl (5) 

where U a i is the mixing matrix between mass and flavor eigenstates, F the flux at the Earth, 
F° the flux at production, and the suffix x represents \x and r, of which flavor neutrinos are 
considered to have the same flux. As \U e i\ 2 ~ 0.7 for the LMA models, the flavor mixing is 
expected to harden the v e spectrum. 



1 Here, we note that there is a wrong statement in Section 2.2 in AST. We wrote there "We can naively deal 
with anti-neutrino oscillation effect as vacuum oscillation, since i> e 's are not affected by the resonance." 
As is illustrated in the text, this is wrong, however, the calculations in AST is including supernova matter 
effect appropriately. 
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However, in the case of inverted mass hierarchy, this situation changes dramatically. Since 
v 3 is the lightest, u e is created as v 3 in the supernova core. In that case, it is well known 
that at a so called resonance point, there occurs a level crossing between v\ and v 3 . At 
this resonance point, complete V\ <-» v 3 conversion occurs when the so called adiabaticity is 
sufficiently small compared to one (it is said that resonance is nonadiabatic), while never 
occurs when it is large (adiabatic resonance). In the lower panel of Fig. ^ we show the 
adiabaticity 7 of each model, which is written by 

= Amf 3 sin 2 26 13 
2E u cos26i3\d\nn e / dr\ re s 

where E v is the neutrino energy and n e is the electron number density; we have assumed 
E v = 20 MeV in the figure. 2 The behavior of the adiabaticity shown in the lower panel 
of Fig. ^ can be qualitatively understood as follows. The supernova profile on which the 
adiabaticity depends is the gradient of logarithmic density, \d\nn e / 'dr\~ l . Then, if we assume 
the simple power-law density profile (n e = nor~ a , no and a are constants), the adiabaticity 
is proportional to the radius. This global behavior can be seen in the lower panel of Fig. C] 
changes in the power-law index a and hq are responsible for the variability of 7 in the small 
scale. 

From Fig. d we can expect that for INV-S, the resonance is always nonadiabatic. On 
the other hand, for INV-L, it is basically adiabatic but at the first resonance point, the 
adiabaticity becomes comparable to one particularly at > 10 sec. The situation is very 
complicated for INV-L model, and there is no way but to calculate numerically for the 
case. From here, however, we estimate the flux in the case of the resonance is completely 
nonadiabatic or completely adiabatic, for simplicity. When the resonance is nonadiabatic, 
the situation is the same as in the case of normal mass hierarchy (because u e at production 
becomes V\ at the stellar surface), and the v e flux we observe is represented by Eq. (JHJ). On 
the other hand, adiabatic resonance forces v e at production to become z/3 when it appears 
from the stellar surface and therefore, the observed z/ e flux is given by 

F> e = \U e3 \ 2 F° e + (l-\U e3 \ 2 )Fl. (7) 

Since \U e3 \ 2 is constrained to be much smaller than 1 from reactor experiment [20], Eq. (J2J) 

2 Note that the definition of 7 [Eq. ©] is available only at the resonance points, although in Fig. Q]we 
have shown it as a function of radius over the entire region in supernovae. 
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indicates that complete conversion takes place between v e and v^ T . 



III. RESULTS 



We calculate the SRN flux using the formula 



dF v r 2 - n , .dNJEL) ,dt 7 

C L ^sn(^) :l u) (l + z)-dz, (8) 



dE v Jo dE' v dz 

where E' u = (l+z)E u , R SN (z) is supernova rate per comoving volume at redshift z, dN u /dE u 
energy spectrum of emitted neutrinos, 2 max the redshift when the gravitational collapses 
began, which we assumed to be 5. As supernova rate, we use the most reasonable model 
to date, which is based on the optical/UV observation of star formation history in the 
universe by Hubble Space Telescope |2l|, and the model was also used in AST as "SN1". 
In this model, the supernova rate exponentially increases with z and peaks at z ~ 1.5 and 



exponentially decreases in further high- ,2 region (see Fig. 1 in Ref. |22(]). Since for the 
optical/UV observation, there is still uncertainty due to dust extinction, we introduce a 
factor / as Rs${z) = fRsm(z) and assume that it is independent of z. (This argument is 
almost the same one as that in Ref. So The case / = 1 corresponds to the local supernova 
rate of .Rsn(O) = 8.5 x 10~ 5 h? yr" 1 Mpc~ 3 , while the observational local supernova rate 
is R SN = (1.2 ± 0.4) x l(T 5 /4 yr" 1 Mpc" 3 tJ, where h 70 = H /70 km s' 1 Mpc" 1 . For a 
while, we set / = 1 fiducial value. 

Figure El shows SRN flux for the various oscillation models explained in Section |H] Flavor 
conversions enhance the average v e energy. Three models NOR-S, NOR-L, and INV-S are 
degenerated and the model INV-L has the hardest energy spectrum as is expected from the 
qualitative discussions in Section Hi Bl [see Eqs. © and (J7J)]. In the same figure, we also 
show the flux in the case of no oscillation, for comparison. In that case, the flavor conversion 
does not take place, then the original z/ e flux is actually observed [F Pe = F® e , in contrast 
with the Eqs. (jSJ) and ((ZJ) for the oscillation cases], resulting in the softest spectrum in the 
five models under consideration. 

n 

Recent SK upper bound for the SRN flux p| constrain the theoretical predictions by 
factor ~ 3. For example, AST predicted the SRN flux integrated over entire energy to be 
11 cm" 2 s -1 , while the corresponding SK limit calculated with the AST spectral shape is 
31 cm -2 s _1 . The AST model cited above corresponds to NOR-L model in this study. We 
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repeat the analysis given in the SK paper [3] and obtain the flux upper limit at 90% C.L. for 
various oscillation models. We summarize the results in Table [I] As shown in the table, all 
the oscillation models are not ruled out yet, since the theoretical predictions are still smaller 
than the corresponding SK limit, while the observational upper limit is rather severer for 
the INV-L model. 

TABLE I: The predicted SRN flux for various oscillation models and the corresponding SK limit 
(90% C.L.) {J. The ratio between the prediction and the limit is shown in the fourth column. 
Event rate from 19.3 MeV to the energy where the SRN flux dominates the atmospheric v e flux is 
also shown in the fifth column. 

Model Predicted flux SK limit (90% C.L.) Prediction/Limit Event rate for E v > 19.3 MeV 



NOR-S 


12 cm" 


2 s" 


-1 


< 35 cm" 


" 2 s" 


-1 


0.34 


0.80 /year (E u < 30 MeV) 


NOR-L 


11 cm" 


" 2 S" 


-1 


< 34 cm" 


" 2 s" 


-1 


0.33 


0.81 /year (E u < 30 MeV) 


INV-S 


11 cm" 


" 2 S" 


-1 


< 34 cm" 


" 2 S" 


-1 


0.33 


0.81 /year (E u < 30 MeV) 


INV-L 


9.0 cm" 


" 2 s- 


-1 


< 12 cm" 


" 2 s" 


-1 


0.74 


2.0 /year (E u < 37 MeV) 


no oscillation 


12 cm" 


" 2 s" 


-1 


< 73 cm" 


" 2 s" 


-1 


0.17 


0.43 /year [E v < 27 MeV) 



IV. DISCUSSION 

A. Comparison with the other calculations 

The conversion probability for INV-L model is expected to be time-dependent as discussed 
in Section III Bl We show in the upper panel of Fig. El the conversion probability P{p e — ► 
v^) obtained with density profiles for 5, 10, and 15 seconds after bounce. The conversion 
probability in the case that the static progenitor model is used is also shown for comparison. 
At 5 seconds after bounce, it is almost the same for all energy as in the case of static model, 
while at 15 seconds, it changes by ~ 10% at 30 MeV, and by ~ 20% at 60 MeV. However, 
since almost all of neutrinos are emitted by 5 seconds after bounce, the effective conversion 
probability is expected to be almost the same as that in the case of static model. Actually, 
we calculated the time-averaged probability by the flux and the result is also shown in the 
same figure, indicating only ~ 5% difference even at 60 MeV. 

We also compared the flux calculated in Section IIHI with that obtained by the AST 
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TABLE II: Flux difference from the model obtained with static progenitor model and without the 
Earth matter effect. The flux is integrated by energy. The second column is for the flux integrated 
over the entire energy range and the third column is over the detection energy range, or E v > 19.3 
MeV. The + (— ) sign of each entry represents that the flux increased (decreased) by the following 
number. 



Model 


Total flux 


Flux for E v > 19.3 MeV 


NOR-S 


+0.015% 


-4.0% 


NOR-L 


+0.015% 


-3.9% 


INV-S 


+0.015% 


-4.0% 


INV-L 


-0.43% 


-2.4% 



calculation |2|. The AST calculation adopted the static progenitor model and did not 
include the Earth matter effect. We show in the lower panel in Fig. El the flux difference 
as a function of neutrino energy. For the three models NOR-S, NOR-L, and INV-S, the 
difference comes from the presence of the Earth matter effect, since these models are not 
affected by the time-dependent density profile as illustrated in Section III Bl The Earth 
matter effect changes flux by ~ 15% at 60 MeV. For INV-L model, the Earth matter effect 
is absent (for this reason, see e.g., Re, 0). and th us, the fluX difference comes from the 
difference of the supernova model (static or time-dependent). As shown in the upper panel 
of Fig. El this effect is small and it makes only ~ 5% difference at 60 MeV. In Table ITTl we 
show the difference of the energy-integrated flux. In the detection energy range, E v > 19.3 
MeV, a few % difference is expected. 



B. Implications for neutrino oscillation and the future 

Theoretical predictions of the SRN flux depend on the factor /, which we have assumed to 
be one until here, in a proportional manner. This parameter / is not known well because it 
concerns dust extinction of optical/UV photons, and actually even though we take / = 2, it 
does not conflict with the local supernova rate observation. [Observational local supernova 
rate is R SN = (1.2 ± 0.4) x 10" 4 tf yr" 1 Mpc" 3 Q, while SN1 model in AST predicts 
-Rsni(O) = 8.5 x 10~ 5 hj yr" 1 Mpc~ 3 .] However, even if we adopt / > 1.4, we cannot 
conclude that the INV-L model is already ruled out, although from Table |l] it appears that 
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the predicted flux of the INV-L model becomes larger than the corresponding upper limit. 
This is because the SRN flux contains uncertainty other than the parameter /, which is that 
concerning the original neutrino spectrum: in particular, the difference of average energies 
between i? e 's and z/ x 's, about which we give a more detailed discussion in the next subsection. 
Thus, although the current SRN observation includes useful information on the combined 
quantity of supernova rate in the universe and neutrino oscillation parameters, it is difficult 
to give an explicit solution to the one specific problem at the present stage. 

Now, we consider the future possibility to detect SRN or to set severer constraint on 
neutrino models. The largest background against SRN detection at SK is so called invisible 
muon decay products. This event is illustrated as follows. The atmospheric neutrinos 
produce muons by interaction with the nucleons (both free and bound) in the fiducial volume. 
If these muons are produced with energies below Cherenkov radiation threshold (kinetic 
energy less than 53 MeV), then they will not be detected ("invisible muons"), but their 
decay-produced electrons and positrons will be. Since the muon decay signal will mimic the 
u e p — > e + n processes in SK, it is difficult to distinguish SRN from these events. Recent 
SK limits are obtained by the analysis including this invisible muon background. In the 
near future, however, it will be plausible to distinguish the invisible muon signals from the 
SRN signals, using the gamma rays emitted from nuclei which interacted with atmospheric 
neutrinos 24|. Therefore, if we can detect gamma ray events, whose energies are about 
5-10 MeV, before invisible muon events by muon life time, we can subtract them from the 
candidates of SRN signals. In that case, the upper limit would be much lower (by factor 
~3) when the current data of 1,496 days are reanalyzed 24|, and the SRN signal might 
be detected or the more powerful information on the combined quantities of the supernova 
rate and the neutrino mixing parameters would be obtained. In Table HI we also show the 
expected event rate of the SRN signal for various models. The integrated energy ranges are 
set from 19.3 MeV to the energy where the SRN flux dominates the atmospheric v e flux. 
Then, without invisible muon events, the SK data of 1,496 days would be sufficient to permit 
the SRN detection and set very severe constraint on the neutrino mixing parameters. In 
particular, inverted mass hierarchy with large # 13 (INV-L model) would be ruled out first 
among five models we have considered. 

There is another possibility to enhance the average v e energy, or resonant spin-flavor 
conversions. This mechanism is induced by interaction between nonzero magnetic moment 
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of neutrinos and supernova magnetic field. We investigate this mechanism in another paper 
2511 for details. 



C. Uncertainty concerning supernova neutrino spectra 



In this study, 



ermore group 



we used the original neutrino spectrum calculated by the Lawrence Liv- 



11 1 . Unfortunately, their study as well as the other published full numerical 
supernova collapse simulations have not yet included the nucleon bremsstrahlung process or 
nucleon recoils, even though it is no longer controversial that these effects are important. 
Recent studies (e.g. Ref. [12j) including all these processes have shown that average u x 
energy exceeds the average i> e energy by only a small amount, 10% being a typical number. 
If it is the case, the oscillated SRN flux is likely to be close to the value in the case for no 
oscillation, and the recent SK limit is not severe enough to constrain the various oscillation 
parameters. However, it is premature to conclude that their results are correct, since it is 
based on the neutrino transport study on the background of an assumed neutron star atmo- 
sphere, and this approach lacks hydro-dynamical self-consistency. Further, it is also because 
the mean energies and their ratios change significantly between the supernova bounce, ac- 
cretion phase, and the later neutron star cooling phase. Whichever is the case, the future 
SK limit is expected to be sufficiently severe to constrain various oscillation parameters as 
well as supernova rate in the universe. 
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Radius | R Q ] 

FIG. 1: Upper panel: The pY e profile of supernova at 0.5, 2, 5, 10, and 15 seconds after bounce. The 
horizontal band shows resonance condition for two INV models (the band width comes from the 
energy range 5-70 MeV), i.e., at intersections between the pY e curve and the horizontal band, the 
MSW resonance occurs. Lower panel: The adiabaticity of the resonance (i.e., if the adiabaticity 
at the resonance point is larger than 1, then the resonance is adiabatic) for INV-L (labeled by 
sin 2 26*13 = 0.04) and INV-S (sin 2 20i 3 = 10" 6 ) models. The neutrino energy is assumed to be 20 
MeV, and the line types are the same as those used in the upper panel. 
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NOR-S 
NOR-L 




Neutrino Energy [MeV] 



FIG. 2: Number flux for zv e 's for various neutrino oscillation models. The spectra for NOR-S, NOR- 
L, and INV-S are degenerated, while that for INV-L is the hardest one. The flux of atmospheric 
v e is also shown. 
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FIG. 3: Upper panel: Conversion probability P(v e — > t^) as a function of neutrino energy calculated 
using static progenitor model (dot-dashed line), shock model at 5 (dotted line), 10 (short-dashed 
line), and 15 seconds (long-dashed line) after bounce. Time-averaged probability by the flux is also 
shown by solid line. Lower panel: Flux difference from the model obtained with static progenitor 
model and without the Earth matter effect, for various oscillation models. 
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